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Abstract 

The list of properties of Fibonacci numbers F„ (with multifaceted relevance 
in physics) is complemented by an empirical observation that in combination 
with the "next" family of the "delayed Fibonacci" numbers Gm (called, for 
convenience, "Gibonacci numbers" here), both sets exhibit certain remarkable 
and fairly unexpected asymptotic mutual-bracketing properties. 
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1 Introduction 



Our forthcoming considerations were inspired by the decoration of the north- 
most subway station "Holesovice" in Prague. Any passenger who waits there 
for the train to the city may wonder why the architects ornamented the walls 
by the regularly repeated hnes of 2, 2, 3, 4, 5, and 6 tiles. In this sense, our 
present paper comes only too late when suggesting a replacement of such a 
pattern by an inessentially modified sextet of the integer numbers = 2, 
Gi = 2, Gs = 3, Ge = 4, = 5 and Gg = 7. 



2 Fibonacci-type recurrences with zero- and 
one-step delay 

Our deeper motivation stems from the fact that the above-mentioned segment 
of a sequence of G„ may be understood as one of the most natural "delayed" 
modifications of the famous Fibonacci numbers [1]. The latter sequence with 
conventional denotation F{n + 1) = F„ [2] is defined by the well known and 
extremely elementary recurrent relations 



Fn — Fn-l + Fn-2, 



n = 2,3. 



Fo = i^i = 1 



(1) 



Similarly, their "delayed" version may be generated by the very similar "next- 
to-Fibonacci" recurrences with a one-step delay. 



Gn — Gn-2 + Gn-z, 



n = 3,4, ... , 



Go = Gi = G2 = 1 . (2) 



In our present note we intend to present a few arguments showing why both 
the sequences F„ and G„ might be considered comparably interesting. 

Firstly, let us remind the reader that both these sequences are virtually 
equally easy to construct. Both difference equations (1) and (2) for F„ and 
Gjn have constant coefficients and may be analyzed by the similar standard 
ansatzs, viz.. 



and 



Gr,^A 



Q{A) 



Via) 



B 



+ b 



Q{B) 



c 



Q{C) 



(3) 



(4) 



respectively. This procedure leads to the similar implicit definitions of the 
necessary quotients. 
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2.1 Fibonacci sequences 



In the former case which corresponds to the usual Fibonacci numbers F„ the 
imphcit definition of the quotient is shghtly simpler, 



V{a,b) = V{a,b) + 1 ■ 



(5) 



This equation may be immediately assigned the following well known and com- 
plete explicit solution, 



Via) 



l + \/5 



1 - ^/5 



2 ' 2 
In it, the former item coincides with the famous "golden mean" value. 



(6) 



2.2 "Gibonacci" sequences 

There is really no reason why one should be afraid of searching for the similar 
explicit solutions in the delayed case with the similar implicit definition 

-|3 

Q{A,B,C)\ = Q{A,B,C) + 1 (7) 

of the quotients. Routinely, one arrives at just one real and positive root 



Q^A) ^D+ + D. 



(8) 



where we abbreviated 



D^^ ('1±^V^' ~ 0.98699 , 



D_ = ('i-^V^'« 0.3377 



and where 



/23 

xo = W — 0.9229582. 
V 27 



The other two complex conjugate roots possess the equally compact represen- 
tation 



Q{B) = [Q{C) 

where (p = 27r/3 and q^d) = D+ - D_ 0.649264. 



e^-D^ + e-^-D_ (= g^^^ + g^^^ (9) 
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2.3 Specific initial boundary conditions 



Once we require the compatibility of formula (3) with the Fibonacci's boundary 
conditions Fq = Fi = 1, we have to extract the corresponding values of a and 
b from the two equations 



1 l + \/5 , ,1-^5 
1 — a + — a — h 



2 "2 

This gives the well known explicit formula for Fibonacci numbers, 



(10) 



(11) 



In the similar manner, the delayed-Fibonacci general solution (4) comple- 
mented by the boundary conditions Gq = Gi = G2 = 1 will define our present 
Gibonacci numbers Gm- First of all we must find the correct values of the 
coefficients in eq. (4) by solving the the triplet of the linear equations 



l^A + B + C^Ag^A)+BQ^B)+Cg^c) = A 



QiA) 



+ B 



n 2 



QiB) 



+ G 



n 2 



(12) 



Its numerical solution gives the approximate form of the result. 



A = 0.72212441830311284114 . . . , 
B = 0.13893779084844357942 ...-i- 0.20225012409895253966 . . . , (13) 
C = 0.13893779084844357942 + 0.20225012409895253966 .... 

Now, the same result will be derived in closed form, completing in this way the 
analogy with the previous "non-delayed" formula (11). 

Firstly, we suppress all temptations to use a computerized symbolic manip- 
ulations and put B — K — i L — C*. The first item in eq. (12) then offers the 
direct elimination of ^4 = 1 — 2K. Moreover, after an abbreviation 

l+2/s:(cos(^-l)+2L sin(^ = E, 1+2 /sT (cos<^-l)-2 L sin 93 = A , (14) 

we may re- write the remaining two lines of eq. (12) in a particularly friendly 
form, 

E + A D_ = 1 



3 E ^2 ^ 3 A = 1 . 



(15) 



This two-by-two matrix equation is readily solvable. 



3x0 
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The final backward insertion in eq. (14) is trivial and gives the final answer, 
6X = 2 - (D+ - /}_) {3D+ + 3L>_ + 1) , 

OXo 

2V3L^ — ( 3Dl + 3Dl - D+ - dJ) 
which is compatible with its numerical check (13). 



3 Intertwining behavior of the Fibonacci-type 

sequences 

3.1 Asymptotics of F„ and Gm at the large indices 

As long as wc observe that |?7(a)| = 1.618... > 1 and \ri{h)\ = 0.618... < 1 
while = 1.3247. . . > 1 and = = | - 0.662 . . . ±^0.562.. . |, 

i.e., \g(^B)\ = .8688... < 1, wc may conclude that irrespectively of the initial 
boundary conditions the solutions of both the difference equations (2) and (1) 
will always exhibit a similar asymptotic behavior. The reason is that both Fn 
and Gm are dominated by the single power-law term, i.e.. 



Fn^a 



+ n>l (16) 



while 



"1 m 

Gm — A g(^A) 



+ m>l. (17) 



Prom these relations we may deduce that the size of the numbers and Gm 
cannot remain comparable unless the indices N ^ 1 and M ^ 1 obey the 
following rule, 

^-^ = 0.584.... (18) 
M lnr7(,) ^ ' 

This means that the replacement of the Fibonacci recurrences (1) by their one- 
step-delayed modification (2) slows down the asymptotic growth of the new 
sequence, Gm ~ -^entier{o.584M} at M > 1. 

3.2 Inequalities between Fn and Gm at the finite indices 



An identity 

= 0.58435715765740408667... = — + A ( — V+ A— , 
ln77(„) 12 103 V12; 106 109 VlOii 

(19) 
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complements eq. (18) and indicates that the ratio 7/12 may play a key role 
in our present analysis. Indeed, once we tentatively re-index all the very large 
Fibonacci numbers, 

Fj,^f(j,k), N = N{j,k) = 7j + k, 0<k<6, (20) 

and once we re-write their delayed-generated alternative in the similar form, 

GM = g{J,K), M = M{J,K) = 12 J + K, < K < 11, (21) 

we may re-interpret the above "asymptotic comparability rule" (18) as a re- 
quirement 

+ ^ - ^ (22) 



12J + K 12' 

In the other words, we achieve the approximative asymptotic coincidence of 
Gm = gj,K with F/v = fj,k if and only if the new auxiliary indices j and J do 
not differ too much. 

3.3 Inequalities between F„ and Gm at the small indices 



A priori, there is no reason to believe that the similar rule could be extended 
to the domain of the small indices. Nevertheless, we may take the union set of 
all the numbers in eqs. (20) and (21) and orc?er this family in a way starting 
at the very first subscripts. In this way the first seven lines of inequalities are 
revealed. 



9(0, -l] 


(=0) </(0,0) ( 


= 1) < 9(0,0) ( 


= 1 


9(0,1) ( 


= 1) </(0,l) (= 


1) < 9(0,2) (= 


1), 


9(0,5) ( 


= 2) </(0,2) (= 


2) <^?(0,4) (= 


2), 


^?(0,5) ( 


= 3) </(0,3) (= 


3) < (7(0,6) (= 


4), 


9(0,Q) ( 


= 4) </(0,4) (= 


5) <^7(0,7) (= 


5), 


9(0,8) ( 


= 7) </(0,5) (= 


8) < (7(0,9) (= 


9), 


^(0,10) 


(= 12) </(0,6) ( 


= 13) <^(0,11) (= 



(23) 
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Encouraged by the smoothness of this pattern we may verify, with an utterly 
unexpected success, the existence and vahdity of its next-step continuation 



^(0, 11) (= 16) < /(I, 0) (= 21) < g{l, 0) (= 21), 



5(1,1) 
5(1,3) 

5(1,5) 

5(1,6) 

5(1,8) 



= 28) </(l,l) (=34) <^(1,2) (=37), 
= 49) </(l,2) (=55) <^(1,4) (=65), 
= 86) </(l,3) (=89) <^(1,6) (=114), 
= 114) < /(1, 4) (= 144) < ^(1,7) (= 151), 
= 200) < /(I, 5) (= 233) < ^(1,9) (= 265), 



(24) 



5(1, 10) (= 351) < /(1, 6) (= 377) < 5(1, 11) (= 465) . 

Now, there comes one of our main empirical observations. Against all odds, 
the same scheme works during unexpectedly many iterations numbered by the 
integer K which appears as the first argument in the functions g{K,*) and 
f{K, *) and which was equal to zero in eq. (23) and to one in the subsequent 
set of the fourteen inequahties (24). We arrive at the formidably extensive set 
of the empirical inequalities 

g{K -1,11) < f\K, 0) < g{K, 0) , < 48 

g{K,l) <f{K,l) <g{K,2) , K < 35 

g{K,3) <f{K,2) <g{K,4) , K < 21 

g{K,5) <f{K,3) <g{K,6) , K<7 (25) 

giK,6) <f{K,A) <giK,7) , K < Al 

g{K,8) <f{K,5) <g{K,9) , < 27 

g{K, 10) <fiK,6) <g{K, 11) , K<U. 

The existence and structure of the upper limits of their validity reminds us of 
the fact that the next correction to cq. (22) (given in eq. (19)) does not vanish. 
This means that the range of the allowed K in eq. (25) cannot be unlimited. 



3.4 Inequalities between Fn and Gm at the growing in- 
dices 

It is remarkable that the inequalities (25) are violated so extremely slowly and 
in such an unbelievably regular manner. This is one of the main consequences 
of the smallness of the absolute value of the second correction 0.003 ■ (7/12)^ 
to the rule N/M — 7/12 because it will enable us to extend our inequality 
pattern beyond its limits listed in eq. (25). 

As long as the next correction to the asymptotic comparability rule N/M — 
7/12 is positive, we know in advance that the middle terms f{K, J) in eq. (25) 
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will all move to the left and violate their lower estimates at a critical K. Due 
to this expectation (confirmed by the explicit calculations in MAPLE), the 
next-step bracketing law acquires the mere following shifted form 



g{K -1,10) <fiK,0) <g{K- 1,11) 



g{K,Q 

9{K,A 

g{K,h 

9{K,1 

g(K,9 



<f{K,l) <g{K,l) 
<f{K,2) <g{K,3) 
<f{K,3) <g{K,5) 
<f{K,A) <g{K,6) 
<f{K,5) <g{K,8) 
<fiK,6) <giK,10) 



48 < X < 96 
35 < X < 82 
21 < X < 70 
7 <K <55 
41 < X < 89 
27 < X < 75 
14 < X < 61. 



(26) 



Etc. One should notice that due to the fact that all the K^s arc ahcady 
very large, all the separate intervals of validity of the innovated rule (26) arc 
perccivably longer than their respective predecessors in eq.(25). This seems to 
indicate a general tendency, the detailed analysis of which would require much 
more space than available here. 



4 Outlook 

4.1 Towards combinatorial applications 

Fibonacci recurrences (1) (without any delay) are extremely popular and find 
(perhaps, unexpectedly) numerous practical applications. Pars pro toto, Fi- 
bonacci numbers F„ occurred recently as a sequence which numbers all the 
possible re-arrangements of the Born-Lanczos expansions of the scattering am- 
plitudes in quantum mechanics [3] . The mathematical essence of this particular 
application lies, in a way illustrated by Table 1, in the following combinatorial 
representation of the Fibonacci numbers. 

One may feel inspired to re-interpret the latter property as a definition. In the 
next step one then could modify this type of definition, obtaining the following 
"higher" Fibonacci numbers, 

'"=„.£.,('-/')■(:>(';■)*(';*)-" 
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(cf. Table 2) or 

(cf. Table 3) etc. From this background, one can always return to the recurrent 
approach, discovering that it starts from a general initial £— plet of values 

= k^O,l,. ..,£-!,£, £=1,2,.... (30) 

Moreover, the old combinatorial definitions (28) or (29) (etc) become replaced 
by the corresponding Fibonnaci-type new three-term-like recurrences 

Fi'^-Fi'}, + Fi%, n = 3,4,... (31) 

or 

Fi'^-Fi'}, + Fi%, n = 4,5,... (32) 

etc. In this way we would obtain a new series of generalized Fibonacci numbers. 
The analysis of the influence of the delay in the underlying recurrences lies 
already beyond the scope of our present short communication but it is worth 
mentioning that it would proceed precisely along the lines applied here to the 
simplest delayed case (2). 

4.2 An appeal of linearity 

An apparent ambiguity of an initialization of the delayed Fibonacci recurrences 
(2) is just fictitious. Indeed, although a more general choice of the initialization 
appears admissible, 

Go{a) = G2{a) = 1, Gi{a) = a e (-oo, oo) (33) 

the question has an elementary answer since the innovated initialization (33) 
merely produces the sequence with elements ^3(0) = 1 + a, ^4(0) = 1 + a, 
G^ia) = 2 + a, ^6(0) = 2 + 2a, ^7(0) = 3 + 2a, Gs{a) = 4 + 3a etc. We 
immediately see that we have 

Gn{a)^Gn-2{'^) + aGn-3{l) (34) 

so that all the variations of a 7^ 1 do not induce any real gain in generality. 

The rule of this type may be understood as one of the manifestations of the 
linearity of our present example. This property opens a path towards appli- 
cations of the similar models in statistics where certain generalized Fibonacci 
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numbers proved related to the close-packed dimers on non-orient able surfaces 
[4], to the local temperature distribution on quasiperiodic chains [5], to the 
existence of fractional statistics in quantum gases of quasiparticles [6], to the 
statistics born by the stacking of squares on a staircase [7] and to the electron 
and phonon excitations in quasicrystals [8]. 

Closely related use of the linearity of the three-term Fibonacci-like recur- 
rences helped, in [9], to introduce randomness directly in the coefficients a in 
a way resembling eq. (34), with possible impact ranging from the modelling 
of chaos (e.g., in quantum gate networks and quantum Turing machines with 
the Turing head controlled by a Fibonacci-like sequence of rotation angles [10]) 
till explicit models of the critical level-spacing distributions, with eigenvectors 
lying between extended and localized [11]. Last but not least, a strong appeal 
of all these models (which may all be interpreted as various discretized ver- 
sions of Schrdinger equations) lies in their generic non-Hermiticity (cf., e.g., 
[12]) which became subject to an intensive study recently (cf., e.g., [13] or all 
papers in the dedicated issue [14]). 

4.3 Numerical "miracles" and open questions 

An exceptional role of our most elementary modification (2) of the most popu- 
lar FihonaccVs recurrences has been illustrated here with a particular emphasis 
on some of the most interesting numerical features of the "Gibonacci" sequence 
Gn- In this sense, our key message has been based, in essence, on the remark- 
ably quick convergence of the asymptotic scries (19). 

We did not pay attention to all the properties of Gm of the similar numerical 
type. For example, we did not throw any light on an alternative numerical 
relation between our asymptotic quotient g(^A) and the golden mean r](^a) which 
is based on the evaluation of the logarithm of their ratio. 

In I ^ I = 0.20001225073664160 .... (35) 

We see that it has an exceptional form (with several zero digits in it) as well as 
a remarkably compact approximate representations with higher precision, e.g., 

- + ^ + — ^ - ^ = 0.20001225073664074 .... (36) 

5 8-104 4-106 27-108 1013 ^ ^ 

The smallncss of the subsequent corrections as well as the use of the natural 
logarithm in eq. (35) do not have in fact any natural explanation. In the other 
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words, the comparatively high reUabiUty of the estimate 

(37) 

represents an unclarified numerical mystery. Why the two quotients ^(a) and 
ri(^a) should be related at all? And even if yes, why are they related just to the 
base e = exp(l) ^ 2.718 of natural logarithms? 

Marginally, we may add a remark that in the next possible study of the 
doubly delayed recurrent relations of the above Fibonacci type, 

Hn — Hn-3 + Hn-A (38) 

one could ask why precisely the one-step delay in eq. (2) should be considered 
exceptional. The first answer could be purely pragmatic, stating that algebraic 
equations (5) and (7) seem to be the only sufficiently easily manageable pair of 
definitions of quotients. Indeed, in the next case with a double-step delay, even 
the straightforward asymptotic analysis could be marred by the less transparent 
solution of the quartic analogue of eqs. (5) and (7), 

= T + 1 . (39) 

Its two real roots 



-0.72449195900051561159 . . . , 1.2207440846057594754 . . . 

and their two complex conjugate partners 

T(±) = -0.24812606280262193189 ...±i 1.0339820609759677567 . . . 

can be hardly expressed via reasonably compact formulae. Moreover, on a 
deeper level one reveals that the absolute value \t[±)\ ~ 1.0633 of the two 
complex roots is bigger than one. This implies that in the doubly delayed case, 
the subdominant component of the asymptotics would not decrease anymore, 
with all the possible related complications which were not encountered just in 
the two above-listed cases, viz., in sequences generated by eqs. (1) and (2). This 
adds a further background to our belief that only the famous Fibonacci numbers 
Fn and their present one-step-delayed "Gibonacci" numbers deserve really 
an exceptional attention. 
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Table 1. Chains of elements 4 with pairwise confluences 'v' (cf. [3]). 



length 


eligible structures 


1 












2 


44 










3 




4^ 


^4 






4 




44^ 


4^4 


^44 




5 




444^ 


44^4 


4^44 


^444 




4^^ 




^(^54 






6 




4444^ 


444^4 


44^44 


4^444 




^4444 


44^^ 


4^4^ 


^44^ 


4^^4 




^4^4 


^^44 
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Table 2. Chains of elements ^ with confluences in triplets X- 



length 


eligible structures 


1 


4 












2 














3 




4 










4 




44 


44 








5 














6 












44 


7 






44444 


44444 


44444 


















8 






444444 


444444 


444444 






444444 


444444 












^^^^^^^^ 
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Table 3. Chains of elements ^ with quadruplets <0'. 



length 


eligible structures 


1 












2 












3 












4 













5 




40 


04 






6 




440 


404 


044 




7 




4440 


4404 


4044 


0444 


8 




44440 


44404 


44044 


40444 




04M4 


00 








9 




444440 


444404 


444044 


440444 






044444 


400 


040 


004 
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